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ABSTRACT 

The generally difficult problem of a system with time 
delay is studied by using parameter space methods and com- 
puter analysis. The solution is proved to be very easily 
obtained and the results presented are extendable under 
certain circumstances to other systems and useful in other 
fields of analysis. 
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INTRODUCTION 



Dead time or transport lag, defined as a fixed time de- 
lay in the transmission of information between two segments 
of a system, is a phenomenon that occurs frequently in every 
field of engineering. The synthesis of a control system 
where dead time tends to occur has most discouraging aspects 
and the first thoughts of the engineer confronted with the 
problem should be how to go about changing the system design 
so as to eliminate its sources. 

However, many times this is not possible due to physical 
limitations in the system. The dead time may be due to hav- 
ing the measuring element and the controller actuator neces- 
sarily separated by a certain distance, as in the extreme 
case of the remote control of a space vehicle. In other 
cases, there is a. delay between the time at which the con- 
troller actuator is moved and the time at which the plant 
reacts to the change, as in the case of an engine-governor 
system. 

Dead time deteriorates seriously the stability of a 
system since it creates a phase lag that constantly increases 
with 03. 

The analytic treatment of the problem is strenuous and 
obtaining an exact solution is very difficult, if not impos- 
sible. The Characteristic Equation (C.E.) of a system with 
time delay has transcendental nature and a Maclaurin series 

s 
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of the exponential term shows that the system has an infinite 



number of roots. 

Bode and Root Locus analysis are practical and efficient 
methods in the study of linear systems, but in the case of a 
system with transport lag they become too cumbersome. They 
are also single parameter techniques. However, it is pre- 
cisely in the case of systems with transport lag that the 
’’Parameter Space Method” (1) , (2) , appears as a most useful 

and practical aid. This method has application whenever the 
system has several adjustable parameters and consists gener- 
ally in the mapping of the S plane in a N dimensional space 
of the parameters of the system. With this mapping completed 
it is possible to choose the N parameter values that allow 
the system to have a certain set of root locations. 

As previously mentioned, a system with transport lag has 
a C.E. with an infinite number of roots; therefore, if we 
were attempting tO' look at stability, the knowledge of the 
locations of some of the roots would not be the answer. 
Methods are available for testing the stability of systems 
with transport lag (2) , (3) . In this work we will concen- 
trate on finding the parameter values that guarantee a near- 
optimal behavior of the system recognizing, however, that 
analysis of stability should be performed after this synthe- 
sis phase is finished. 

The study of the general problem was followed by appli- 
cation to a practical example. 
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Consider the 
fer functions of 



II. THEORETICAL BACKGROUND 

two block diagrams of Figure 1. 
the two systems represented are: 



The trans- 



T.F 

T.F 



We see that the 
depend on the locati 
C.E. is the same ind 
in the forward or ba 
of the following stu 
dead time occurs. 

The term exp(-TS 
nature. In fact, th 

-TS 

e 



^ G(S)e~'^^ 

^ 1+G(S) .H(S) .e"^^ 

^ G(S) 

^ 1+G(S) .H(S) .e“^^ 

form of the C.E. of the system does 

on of the dead time block, that is, 

ependently of having the dead time 

ckward path. This allows the exten 

dy to a large class of problems whe 



( 1 ) 

( 2 ) 

not 

the 

block 

sion 

re 



gives to 


this 


equation 


a transcendental 


Maclaur in 


exp 


ansion of 


exp (-TS) is: 


1 - ST + 


ST^ 


3 

- ST + . , 


(3) 



We see then that a system with transport lag has an in- 
finite number of roots which, other than making complete sta- 
bility analysis very difficult, has the effect of making some 
graphical methods, as Root Locus and Bode Diagram, very in- 
appropriate . 

The stability of systems where dead time appears is al- 
ways a factor to consider since stability is seriously 

-TS 

affected by it. Consider that if we examine e through 
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sinusoidal analysis and let s 



= j03, the transport lag term 



becomes 



F(jco) = 



= 1 e 



-j03T 



(4) 



We see then that the transport term has a magnitude of 
1 and a lag phase angle of 03T . The phase lag added by the 
transport lag term increases as O) increases. The system 
may, therefore, be very difficult to stabilize. 

However, we assumed initially that the nature of the 
system is such that all the roots lie in the left half of 
the S plane. The system is then intrinsically stable and 
the problem is to find the values of the adjustable para- 
meters that make the system behave in an optimal manner in 
relation to an established criterion. We concentrated on 
finding the parameter values that guarantee the existence 
of a set of dominant roots for the system, that is, a set 
of roots close to the origin. Stability analysis, using for 
instance, the methods in references (2) and (3), should be 
used after the parameters for the system were selected. To 
determine the parameter values, the parameter space method 
can be used with great -advantage since, as will be shown, 
it allows us to relate precisely parameter values with root 
locations. We considered the case of systems with three 
adjustable parameters . 
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lag 



The general form of the C.E. of a system with transport 



Is : 



N 



N 



F(S) 




X 



z 




j 



(5) 



For 



i=0 

S = —0 + j (0 



e 




This transformation allows the separation of the C.E. 
into two equations, one containing the real and the other 
the imaginary part. That is: 

F(a,w) = F^(a,03) + jF^(0,o)) = 0 
implies that : 



If we consider that the system has three adjustable parameters 
A, B and C, then they appear in F^(a,03), and/or F^(0,(jo). 

Fixing one of the parameters (for instance, parameter A), 
equations (7) and (8), may be solved for the two parameters 
B and C. Recapitulating, we have that equations (7) and (8) 
are functions of the variable parameters A, B and C of 0 and 
(0 from the assumed roots -0 + jo). Fixing all these variables 
but two of the parameters, for instance, B and C, simulta- 
neous solution of the two equations allows the determination 
of the two values B and C , that together with the fixed 
value of A^ guarantees that the system has a pair of conju- 
gate roots at -0 + jco . 



F (0,(o) = 0 
e 

F (0,w) =.0 
o ’ 



(7) 



( 8 ) 



o 



o 
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Repeating this process for all different values of O) in 

a certain range permits the determination of different set- 

pairs, as in Figure 2* Of course, there is an infinite 

number of these solutions. The union of these solutions 

B.C is a line in the parameter plane of axis B and C as 
o o r r 

is shown in the Figure 3 (online). Working similarly, but 
with the variation of O instead of 0), it is possible to 
obtain similar lines in the parameter plane B-C for which co 
is fixed (o) lines). 

The intersection of a certain line of fixed o) (ca line) 

o o 

and another line of fixed O (o line) gives a point in the 

o o ^ ^ 

parameter plane from which it is possible to read the pair 

ic ic 

of parameters B^ that guarantee the existence of a pair 

of roots for the system of the form -0^ + Figure 4. 

Repeating the above process for different values of A 
allows the determination of a line (or a group of points) 
in the three-dimensional space with ordinates A, B and C for 
which every point corresponds to assigning the system the 
pair of roots -a + io) , See Figure 5. 

If we assume that the system has a third dominant root 
-6^- of a real nature, we can, fixing the value of A, deter- 
mine a line of fixed <5(6^ line) in the parameter plane B-C 
as in Figure 6. This is so because if we substitute the 
solution 6^ in the characteristic equation and fix its value 
together with the value of A, then this equation becomes the 
equation of a line in the plane B-C. 
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Assuming different values for A, but keeping 6 constant, 
the repetition of this process for each value of A results 
in a group of lines in which <S is constant, and these lines, 
if this repetition were performed for an infinite number of 
values of A, form a surface (surface 6^) in the parameter 
space A, B and C for which the system has a real root 6^, 
as in Figure 7. 

In the three-dimensional parameter space, the pierce 
point(s) of the line -cr^ + with the plane -6^ gives one 

or more points that correspond to the set(s) of parameters 



B and 
o 


"o’ 


for which 


the system has a 


set 


of three 


roots 


± 30 ^^ 


and 6 


that we 
o 


assumed would be 


the 


dominant 


ones . 


Notice 


that 


, since we 


are mostly interest 


ed in find 


ing 



roots that lie close to the origin, the range of variation 
of a, 03, and 6 during the process of mapping the S-plane in 
the parameter space does not need to be very large. 

Of course, such an extensive manipulation of equations 
and plot of curves requires the use of a computer. 

We will proceed now to apply these theoretical considera- 
tions to the study of a practical problem, specifically, to 
the determination of the set of near-optimal parameters to be 
used in the operation of an engine-governor system. This 
problem was suggested by the reading of reference (4) and 
subsequent exchange of correspondence with one of its authors 
In this paper, it is suggested that the determination of 
the parameter settings must be made by assuming that the 
dominant roots of the system occupy a. certain well-defined 
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position, specifically that the real and imaginary parts of 
the complex roots must have equal magnitude equal also to 
the magnitude of the third (real) root. 

The parameter space method was used to find new combina- 
tions of parameter values that would permit us to determine 
if the root locations, as assumed in the aforementioned 
paper, would correspond to an optimal solution, and to 
search for new solutions. In addition, the results permit 
study of related problems such as the problem of sensitivity. 
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Ill 



PARAMETER PLANE STUDY OF AN ENGINE-GOVERNOR SYSTEM 



A. GENERAL CONSIDERATIONS 



Most gas engines 


exhibit 


a certain 


degree 


of 


dead 


time 


that is, an interval 


of time 


between a 


change 


in 


fuel 


flow 


and the corresponding 


change 


in driving 


torque 


• 


This 


time 



lag is made of a number of components of which the most 
important are: the delay in "sensing** the speed change that 

follows the load change; the delay within the governing 
mechanism during which the reaction of the sensing element 
is being translated into governor action; the time required 
for the governor action to correct the fuel injection pump 
setting; the delay that must take place before the corrected 
fuel change can be introduced into an engine cylinder; the 
time required to convert the corrected fuel change into 
engine torque. 

Dead time is found to be detrimental to stable operation 
and tends to increase speed transient errors. It also deter- 
mines the maximum rate at which the engine-governor combina- 
tion can be operated. The purpose of this study is to find 
the best isochronous governor parameters so that optimum 
transients, for the operation of gas engines of several 
types be obtained. 

The study of ideal governor settings in the presence of 
dead time will be made in the parameter space. This method 
leads to a computer process for analyzing the variable system 
parameters. The process is restricted to a two parameter space. 
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B. ENGINE ASSUMPTIONS 

The driving torque incremental change AQ is related to 
the small variations of the governor actuator AZ by the 
equation : 

— TC 

(9) 



AQ = Cg e AZ 



— T S 

where e is the dead time operator. 



The minimum 


amount 


of 


dead time T 


that can be 


expected 


given by the 


f ollowing 


formulae : 






T = 


20/N 


for 


a Diesel en 


gine 


(10) 


T = 


40/N 


for 


a two-cycle 


gas engine 


(G.E.) (11) 


T = 


60/N 


for 


a four-cycl 


e gas engine 


(G.E.)(12) 



in which N is the rotational speed in RPM . 

The assumptions made relative to the engine rotating 
parts are that the accelerating torque is equal to the dif- 
ference between the load torque and the driving torque as 
expressed by the following equation: 

TTI 



30 



S AN = AQ - AL 



(13) 



Combination of (9) and (13) gives: 



^ S AN = C, e AZ - AL 
iU o 



(14) 



Change to a system of per unit variables gives the 
engine equation used in the analysis: 

T„ Sn = e”^^Z - Si (15) 

C. GOVERNOR ASSUMPTIONS 

It is assumed that the isochronous pressure compensated 
type of governor is used. The differential equation 
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(16) 



describing its dynamics is: 

(s^ + aj^s)z = 

where “ (N^/Z^)a 2 . 



-p2(S + a2>n 



D. ENGINE GOVERNOR COMBINATION 

Combination of equations (15) and (16) gives; 

(S^+ a^S^+ ^ e"^^S + e"^^)n = - |-S(S+a )il 

M M M 



(17) 



E. USE OF THE PARAMETER SPACE METHOD 

We start with the engine-governor equation (Equation 17) 



(S^+a S^+ ^e“’^^S+ ^Y~^e '‘^^)n = - ^S(S+a )il 
M M M 

From equation (17), the C.E. can be extracted: 



„3 . „2 , ^2 -TS„ . ^2“3 -TS 

S + a-S + — e S + — e = 0 

M M 



(18) 



Since we have three adjustable parameters in the charac- 
teristic equation, we can control the position of three of 
its roots. We assumed then that the system had three domi- 
nant roots similarly to what would happen if the delay in 
the same equation were equal to zero. For these three domi- 
nant roots we had two choices with respect to their intended 
location, that is, we could choose to have three real roots, 
or two complex conjugate roots and a real one. Because we 
were interested in having a small rise time for the tran- 
sients of the system, we chose the second case, that is, we 
assumed that the system should exhibit a pair of complex 
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conjugate roots and a single real root close to the origin 
in the S plane. 

In accordance with this choice we made: 

S = -a + jo) 

We have then that: 

=> (a^ - oj^) - j (2aw) (19) 

= (3ou)^ - a^) + jOa^oj - 03^) (20) 

“TS 

The exponential term e becomes: 

-TS -T(-a+j(i)) 

e = e 

TO “jo)T 
= e e 

To 

= e (cos uT - j sin ojT) (21) 

Substituting (19), (20) and (21) in (18) gives: 



2 3 .. 2 3 ,.. ,2 2 . 

(3ao) - a )+ j (3a o) - u )+ aj^(a - oj - j2aa)) 

^2 Ta 

+ — e (cos 0 )T-j sin uT)(-a+jO)) 

M 

. ^2^3 Ta 



e (cos wT-j sin coT) = 0 



( 22 ) 



M 



Both the real and the imaginary parts of this equation 
must be equal to zero: 

^ ^ T* rv 

3aco^ -a ^+a, (a^-0)^ ) + (o) sin ojT-a cosuT+a^cosojT) = 0 (23) 

M 



2 , ^2 Ta 

3a w-w -2a.au+ ;p~e (a sinwT+co coscoT-a„ s inuT) = 0 



(24) 



Equations (23) and (24) form a system that can give solu- 
tions for any two parameters if the others are fixed. If we 
take and IJ 2 parameters of the plane, and fix the value 

of a^, we can get a complete representation of the loci of 
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the 


points 


(w) , ^ 


■ constant or a^(o), ^2 


(a) , 


0) = 


constant^ 


for which the 


C.E. has a solution of 


the form 


-a + 


jw. 








The curves 


determined as 


above-curves 0 and O) 


are loci 



of all parameter pairs (a^ and ^2) for which the solutions 
-a + jO) and -a-jo) exist. 

If we denote by 6 the dominant real root of the C,E., we 
can by the same methods get a plot of 6 lines in the , VI2 
parameter plane. 

We can repeat the process of finding a, 0) and 6 curves 

for as many values of as we want. 

For a certain = constant, the intersection of selected 

a and 0) curves is the loci of values of the parameters 

and lio for which the selected solution -a + jw exists, 
z o — o 

In the three dimensional orthogonal space and 

a^) the loci of pairs a^. and 1I2 a line along which it is 

guaranteed that any readings of ^3 guarantee the 

existence of the above form of solution. In a similar way, 

the set of lines of constant 6 . one for each value of 

forms a plane in the three 

The single or multiple 

line with the plane of con 

parameters ^2 ^3 ^ 

root 6 and a pair of coni 
o 

since we are primarily int 
set of parameters for the 
meters that offers a *’bett 



dimensional space. ' 

intersections of the aforementioned 
stant ^ gives the set(s) of 

or which the system has a real 
ugate roots -a ^ ^ * However, 

crested in finding the near-optimal 
governor, that is, the set of para- 
er** transient, we will not be 
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looking for a set of the three parameter space parameters as 
indicated above. Instead, we will fix the location of the 
pair of complex roots, which, as explained above, will allow 
the determination of a line in the three dimensional para- 
meter space. We will then proceed to test the system with 
parameters corresponding to points of this line. Of course, 
the parameters corresponding to each of these points deter- 
mine a different location for the dominant real root. 

A computer program was written to plot the a, 03 and 6 
curves on a parameter plane with axis and ]i 2 ^ number 

of values of the parameter (computer program //I) . This 

program was designed with the purpose of making it adaptable 
to studies of any type of engine-governor combination that 
may be described by equation (17). 

We performed the computer analysis first for a hypothe- 
tical system in which the transport lag was zero and later 
for each of the three practical cases cons ider ed--Diesel 
engine, two-cycle gas engine, and four-cycle gas engine-- 
each of them with the same C.E., but with different dead 
times as in equations (10) , (11) and (12) . The problem of 

obtaining parameter values that guarantee certain root loca- 
tions is thus solved. We still needed, however, any kind of 
way of testing the systems with each of the sets of para- 
meters obtained through the parameter space computer program. 
We, naturally, made use of a second computer program (compu- 
ter program //2) which was designed under the following heading. 
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IV. TESTING OF THE ENGINE-GOVERNOR SYSTEM 



We want a program to test each of the systems with dif- 
ferent groups of parameter settings. From equation (17) we 
determine the transfer function (T.F.): 



n_ 



- - s(s + a.) 

M 



o3^ ^2^ ^2 -TS^ ^2^3 -TS 

S + a, S + — e S + — = e 

M M 



(25) 



The denominator of the above T.F. is rearranged as 
follows : 



S3 H. ^ 

M M 

0 ^9 -Tq 

= S^(S + a ) + rr- e (S + a ) 

M 

— TS 

p e (S + a-) 
= S^(S + a^)(l + ^ -2 — 



M S (S + ttj^) 



Then, the T.F. becomes: 



A 



- i S(S + a.) 
M 



-TS 

p e (S + a-) 
S (S +a^) a + — ~2 



MS (S + a^) 



M 



1 + 



-TS.^ ^ N 1+GH 
e (S + a^) 



M 



S (S + a^). 



(26) 



This implies that 



G = - 



TmS 



(27) 



H = - 



(S+a^) 

sTs+op 



(28) 
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Using block diagram description, we can partition G and 



H in the following form: 





And the complete block diagram of the feedback system becomes 
as represented in Figure 8. This configuration allows the 
use of C.S.M.P. 
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V. 



PROCESSING OF A PRACTICAL PROBLEM 



AND ANALYSIS OF RESULTS 

The design of governors requires the choice of a "stan- 
dard" type of response to load changes. Since a change in 
load unavoidably causes a transient change in speed, the 
usual application has the following requirements: 

a) Recovery should be fast and with negligible differ- 
ence in steady state speed; 

b) The maximum deviation from steady state speed should 
be less than some specified percent; 

c) The transient recovery should not overshoot the 
steady state value. 

Trade-offs between these requirements may be necessary in 
some applications . 

For a system with no time delay, but with third order 
characteristics, the best solution for items a) and c) is 
a set of two complex and one real root, such that the real 
root is numerically equal to the real parts of the complex 
roots, and such that the imaginary and real parts of the 
complex roots are equal. For typical applications, a numeri- 
cal choice which gives a satisfactory solution is a real 
root at -3 and complex roots at -3 + j3. (A study of re- 
sults based on other choices is made later.) Justification 
of this is given by Figures 9, 10, 11 and 12, which show the 

transient responses for a system with no time delay, with 
complex roots at -3 + J3 and with several locations of the 
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real root. Notice that Figure 11 is the one that corresponds 
to having the real root 6 equal to -3 and that Figures 9 and 
10 don’t satisfy item a) while Figure 12 shows an overshoot 
which goes against the requirement posed in item c) • 

We wanted to find the near-optimal parameter settings 
for the governor described in reference 4 for each of the 
types considered in the same paper. We assume in all cases 
that the steady state speed of the engine is 450 RPM. 

We used some practical data as it appears in the same 
reference: 

Tm = 5.88 

T = KR/N = 20/450 for a Diesel engine 
= 40/450 for a 2 cycle G.E. 

= 60/450 for a 4 cycle G.E. 

We needed also to know what was the behavior of a system 
governed by equation (17), but with zero delay. This could 
be helpful in establishing the effect of the time delay and 
was used in determining parameter values that, entered as 
data in computer program #2, allowed the plotting of Figures 
9, 10, 11 and 12. 

In each case the parameter plane program was run for 
values of the parameter equal to 0.5, 0.6, 0.7, 0.8, 0.9, 
1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.55, 1.6 and 1.7. 

We decided in all cases that this system should have a 
pair of roots at -3 + j3 and so we concentrated on finding 
the location of the point of the parameter planes where the 
curves co = 3.0 and a = 3.0 intersected. Had the choice of 
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root locations been different, we could still have used the 
same parameter plane plots since they contain curves for 
values of O and o) equal to 0.5, 1*0, 1.5, 2.0, 2.5, 3.0, 3.5, 
4.0, 4.5, 5.0, 5.5 and 6.0. The reason why we did not run 
the program for values of higher than 1.7 was that we 
found that in each case for values of higher than a well- 
determined value, the a = 3.0 and U) = 3.0 curves on which we 
concentrated, would not cross in the parameter plane. This 
is justified by the following: 

For example, in the case of the four cycle gas engine 
with a = 0 ) = 3.0, we have from equation (23) that: 

0 4 

54 + P 2 ^ 0.17 X e * (3 sin 0.4-3 cos 0.4 + a^cos 0.4)= 0 

.254 P2 (0.92 - 1.59)=-54 

- 212.6 

^2 -92a^ - 1.59 

Assuming that the physical nature of the parameters re- 
quires that they have positive values has the effect that 
IJ 2 niust be positive, which implies that 

. 1.59 - - - ^ 

“3 .92 “ ^-72 

The analysis of the parameter plane plots proceeded as 
follows: for each =" const., we determined the intersection 

of the 0 = 3.0 and (0 = 3.0 curves and read the corresponding 
coordinate values of and P 2 * 

Plotting the point correspondent to these parameter 
values in page 2 of the parameter plane plot, (page of curves 
of constant 6) allowed the determination of the location of 
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This value of 



the real root for the considered value of a^. 

6 was further checked by using computer program //3. 

Tables 1, 2, 3 and 4 list in each row the parameter 
values obtained by this process, and that correspondingly 
must be used in testing the dynamics of each of the systems 
so that they present a pair of complex roots at -3 + j3, 
and the value of the real root determined as described. 

.These tables were used to plot Figures 13, 14 and 15 
that represent the projections in each of the orthogonal 
planes that form the parameter space of the points of this 
space that guarantee that each of the systems will have a 



pair o 


f roots at -3 +j3. 


Notice 


in 


these 


pi 


ots that 


the 


points 


correspondent to a 


certain 


sy s 


stem are 


disposed 


in a 


line”. 


Not 


:ice also that 


due to th 


is 


form 


of 


disposi t i 


on , 


if we 


were 


to work with a 


system d 


escribed 


by 


equation 


(17) , 



but with a delay different from those considered here, the 
determination of parameter values that would give the same 
root locations would be possible by simply using interpola- 
t ion . 

Figures 16 through 55 are some of the parameter plane 
plots obtained by using the parameter plane program in the 

t 

IBM 360/65 digital computer. These are only a part of the 
large number of parameter plane plots that we had to obtain 
to construct Tables 1, 2, 3 and 4. We chose those that 
either give parameter values close to those found to the 
"best” or the ones that give real "bad" parameter sets. 
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“3 


“1 


^2 


6 


0.5 


6.6 


126.9 


- 0.60 


0.6 


6.7 


132.3 


- 0.75 


0.7 


6.9 


138.6 


- 0.91 


0.8 


7.1 


145.8 


- 1.09 


0 

VD 


7.2 


150.0 


- 1.29 


1.0 


7.4 


158.0 


- 1.49 


1.1 


7.5 


167.0 


- 1.68 


1.2 


7.9 


177.0 


- 1.95 


1.3 


8.3 


188.0 


- 2 .28 


1.4 


8.5 


198.2 


- 2 . 54 


1.49 


8.97 


209 . 5 


- 3.00 


1.5 


9 . 0 


210.6 


- 3.03 


1.55 


9 . 2 


216.9 


- 3.26 


1.6 


9.3 


225 . 0 


- 3.33 


1.7 


9 . 8 


242 . 0 


- 3.82 



Table 1 - T = 0 
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“3 


“1 


^2 


6 


0.5 


13.4 


193.0 


- 0 .67 


0.6 


12.7 


207.0 


- 0.90 


0.7 


13.4 


225 . 0 


- 1.02 


0.8 


15.4 


251.0 


- 1.25 


0.9 


16 . 8 


280.0 


- 1.58 


1.0 


18.4 


317.0 


- 1.91 


1.1 


21.0 


368.0 


- 2.15 


1.2 


24.1 


438.0 


- 2.54 


1.3 


28.5 


532 . 0 


- 3.00 


1.4 


36 . 8 


707 .0 


- 3 . 69 


1.5 


48.4 


971.0 


- 4.50 


1.55 


65.3 


1300.0 


- 4 . 95 


1 . 6 


89.5 


1859.0 


- 5.50 


1.7 


421.0 


9092 .0 


- 7.20 



Table 4 - T = 60/450 
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Each of the 3 systems was simulated including their res- 
pective time delays and the parameter values as in each row 
of tables 1-4. Figures 56 through 70 are plottings of the 
transients of each of the systems when subjected to 100% of 
load change. We obtained the transient responses for every 
set of parameters in those tables, but we are including only 
as many as needed to illustrate the following analysis. 

Figure 71 is a plot of the transients of each of the 
systems for the case in which the dominant real root has a 
value of -3.0. We can see that all these transient res- 
ponses are very similar. Thus the three roots considered 
in the cases of the systems in which the delay is non-zero, 
are really dominant over all the others, and the determina - 
tion of the near-optimal parameters for a system with trans - 
port lag consists of finding the parameter values that 
determine the location of the dominant roots to be the 
ideal location of the roots for a system with the same 
equation, but with zero time delay . (This criteria is 
restricted to a certain portion of the S plane as will be 
seen later.) 

Examination of Figure 71 shows that for each of the 
systems, the settling time is very approximately the same; 
the differences in time delay determine only the amplitude 
of the transient, this being proportional to the amplitude 
of the delay, which makes sense since the greater the delay, 
the more time the governor takes to interpret the change 
and set the needed corrections. 
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A comparative analysis of Figures 56 through 70 shows 
that the transients of each system with parameter settings 
that correspond to having the real root closer to the origin, 
have larger settling times, while in the opposite case, the 
recoveries are faster, but an overshoot appears. 

The study of the method of determination of parameter 

sets 
of da 
other 
of lo 
F 

that 
part 
t ude 
trans 
words 

sh oul 
excep 

with time delay, this should not apply, and there should be 
some numerical limit on the pattern, depending on the value 
of T . In other words, if we tried to make this root pattern 
fall far to the left of the S plane, either the parameters 
would become negat ive-ther ef or e physically unr ealizable-or 
some other roots would become dominant, or would even fall 
in the right half of the S plane. The answer to this should 



is completed. We tried to make use of the large amount 
ta available by proceeding with the analysis of some 
subjects. The first of these relates to the choice 
cation of the pair of complex roots. 

^ linear system, we know that the same root pattern, 
is, having the real root numerically equal to the real 
of the complex roots and this real part equal in magni- 
to the imaginary part, should give the same type of 
ient, except of course for time scaling. In other 
, the three groups of dominant roots. 




d all give identical transients for a linear system, 
t for the time scale. However, in the case of a system 
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^ m 



1 




be looked at, because, given that this pattern were realiz- 
able more to the left of the S plane, certainly the tran- 
sient would become faster, and the magnitude of the maximum 
deviation could also possibly be smaller. 

We started by finding the sets of parameters that would 
give to the Diesel governor system root sets respectively 
at (-4 and -4 + j4), and (-5 and -5 + j5). The determina- 
tion of parameter sets went as before, that is, in each case, 
we looked in page one of each of the plots of const, for 

the intersections of the 0 = 4 with O) = 4 and 0 = 5 with 
0) = 5 • For each of these intersections, we determined in 
page two of the plot the corresponding value of the real 
root 6, and were finally able to determine the values of 
the parameters and that guarantee the desired 

root locations. The parameter sets determined are displayed 
in the following table: 



DOMINANT 

ROOTS 




“3 


^2 


-4. -4+j4 


16 . 16 


1.84 


444.5 


-5, -5+j5 


20.91 


2 . 25 


643,0 



We used computer program #2 to determine the responses 
of the system while working with these sets of parameters, 
and we obtained the transients displayed in Figures 72 and 
73, Looking at Figure 72, it can be seen that, with the 
parameters set such that the system has dominant roots at 
-4 and -4 + j4, the speed transient is faster and of smaller 



35 



amplitude than any one of those previously found, vhich 



guarantees that the correspondent choice of parameters is 
closer to the "optimal". However, the response, displayed 
in Figure 73 of the system with dominant roots at -5 and 
-5 + j5^ is not good under item c) of the established cri- 
teria. Here the recovery is very fast and the amplitude 
of the transient is even smaller but an overshoot shows. 

This proves that our hypothesis about the existence of a 
limit in the process of bringing the root pattern more to 
the left of the S plane, was reasonable. 

Another topic of interest was the one related to toler- 
ances in the parameter values and/or sensitivity of the 
systems to parameter change. The possibility of and ease 
in performing such studies is one of the great advantages 
of the parameter space method, since the parameter plane 
plots display sets of curves that allow to relate changes 
in parameter values with changes in root location. Figures 
74 and 75 show a parameter plane plot where we marked the 
point that confers to the engine-governor system a set of 
roots at -A and -A + jA. The lines correspond to varia- 

tions of a certain percent in the determined values of the 
parameters and ^2 • possible changes in root locations, 

due to changes in parameter values, is quite obvious. For 
example, the considered variation in the value of while 

]i 2 is kept constant, can cause variation of U) between 1) and 
E and variations in the values of O between B and C. 

Obviously, similar considerations apply to page 2 of each 
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VI 



CONCLUSIONS 



The association of the parameter plane method, with 
appropriate computer analysis, converted a rather difficult 
problem into an almost elementary one. In effect, after 
passing the difficulties in writing and debugging the para- 
meter plane program, the analysis of the information avail- 
able is so easy, that an ill-advised observer may not notice 
all the trouble involved in the solution of the same problem, 
if it were attempted by other methods. 

The problem of finding the set of near optimal parameters 
with time delay consists of finding the parameter values 
that guarantee that the dominant roots of the C.E. of the 
system will be located according to a special pattern. The 
quality of the response improves with the displacement of 
the dominant roots to the left of the S plane but, this pro- 
cess has a limit depending on the value of the delay. The 
determination of this limit at which the quality of the res- 
ponse should be near-optimal has to be made through a trial 
and error process. 

The method described in this study is very easy to follow, 
and the information available allows for investigation in 
many areas of interest, such as sensitivity of the systems to 
variations in parameter values, analysis of permitted toler- 
ances in these values, optimization of performance under a 
defined cost function, and others. 
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Figure 1, Control systems with transport lag in the 
forward and backward paths. 
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Figure 2, Plot of simultaneous solution of equations 
(7) and (8) , _ _ ^ 




Figure 3. 0^ line in the A = const, plane. 




Figure 4. a and u lines in A = const, plane, 
o o 
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Figure 5, 



a . 03 line in A, B, C parameter space, 
o o 




Figure 6. 



a line in A 
o 



const . plane . 




Figure 7 . 



6 plane in A, B, 
o 



C parameter space. 
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diagram of the engine-governor system. 




Figure 9 , Transient response of a hypothetical system 
with zero time delay. 
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Figure lO, Transient response of a hypothetical system 
with zero time delay. 
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Figure 11. 
system with 



"Best” possible response 
zero time delay. 



of 



a hypothetical 
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Figure 12. Transient response of a hypothetical system 
with zero time delay. 
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Figure 13, Projection in the parameter plane of 

points of the parameter space that guarantee that 

each of the systems will have a pair of roots at -3 + j3 
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Figure 15, Projection in the V^2 P^^^^^ter plane of the 

points of the parameter space that guarantee that 

each of the systems will have a pair of roots at -3 + j3. 



49 



•1 






mm 



O'jIi OZ7 Qjh r)S«i Q^3 



zo.s p- 

BtI 

■r 






Pl-o2. 

t 

f:! 



2.p3. 

5^^ 



P3.^. 

si^ 



Rs I 
1 



1 1 
( j 




50 



3'j'} fiir QZ7 Qjfi 0G<1 




51 




52 



0Z7 Q3fi Oib* nua Q53 




53 



Ojj Gj.3 D2V’ 03^ QS3 




54 



riC 02/ 03F Q-i:i QGia qG3 




55 



t»3Q 2^Q UcO ^2C CTQ COO 




56 



027 G3H QG-i 053 072 




57 




I 




CIS 02r 03fi O'^S OGQ Q63 




58 



m 




QO J nxn QZ7 03b 0 OS'! 053 




59 



ESQ bDQ St*C 9EQ 27Q SIC CCC 




60 



0Z7 Q3h UVo OGa Q53 




61 



oa'3 012 03F) oaG ogq qb3 




62 



r:iC 0Z7 053 




63 



C27 Q3fr QVo aS3 




65 



DD3 nic UZT DG<1 qS3 




66 








f 




ao;j 012 027 03^^ a*ro osi qs3 




67 



CIS a^r 03B oat oca qca 




68 



Qj3 OiS 227 23b 2 73 021 




69 



027 03r, QVJ Q53 




70 




71 



Plot of 0, 0) lines. 

T= iO/iSO A3= /J 

X SCALE=06.0 UNITS/INCH 
Y SCALE=90.0 UNITS/INCH 



r)b3 




73 



027 030 Qib' 0G<1 DG3 07Z 




74 




75 



Plot of a, to lines. 
T= ^0/4S0 A3=/.y 

X SCALE=06.0 UNITS/INCH 
Y SCALE=90.0 UNITS/INCH 







: =l|: 



Plot of 6 lines . 



T= 4o/4^o 



AZ= i."] 



X SCALE=06.0 UNITS/INCH 
Y SCALE=90.0 UNITS/INCH 



010 m.7 OGa qea 




78 



rjic Q27 03^ 3 Q53 




79 



QQ*3 C.113 QZT 03P D'^S Qf53 




80 



CilC 027 03h OVi 051 Q53 




81 



CIS DZ7 D3H Diti OSl Q53 




83 



OC’2 CIS 027* 03fi OflG; 05^ Q63 




84 



Q27 031^ a;s 05^1 Q53 




85 



JOii Olfl G27 G3B 05<1 Q63 




86 




I 






DZ7 02^- Qiti 031 Q53 




87 



cm 210 0Z7* D36 Dsa qs3 



.1 







Plot of a, CD lines. 

T=^o/4fo A3=/7 

X SCALE=06.0 UNITS/INCH 
Y SCALE=90.0 UNITS/INCH 



Fig . 



:q 



024 






88 



2Gj UlC QZ7 Q3h Q’iS QS^ QS3 




89 




Figure 56. Speed transient of Diesel. 
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Figure 57. Speed transient of Diesel. 
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Figure 58. Speed transient of Diesel. 
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Figure 59. Speed transient of Diesel. 
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Figure 60. Speed transient of Diesel. 
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Figure 61. Speed transient of 2-cylinder G.E. 
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Figure 62. Speed transient of 2-cylinder G.E. 



96 




Figure 63. Speed transient of 2-cylinder G.E. 
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Figure 64. Speed transient of 2-cylinder G.E. 
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gure 65. Speed transient of 2-cylinder G.E. 




Figure 66. Speed transient of 4-cylinder G.E. 
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Speed transient of 4-cylinder G.E. 




Figure 68. Speed transient of 4-cylinder G.E. 
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Figure 69. 



Speed transient of 4-cylinder G.E. 
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Figure 70. Speed transient of 4-cylinder G.E. 
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Figure 71. Speed transient of the four systems when the 
three dominant roots are located at -3 and -3 + j3. 
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Figure 72, Speed transient of Diesel (Dominant roots 
of system at -4 and -4 + j4). 
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Figure 73. Speed transient of Diesel (Dominant roots 
of system at -5 and -5 + j5). 
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